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The exact diagonalization technique is used to study many-particle properties of interacting elec- 
trons with spin, confined in a two-dimensional harmonic potential. The single-particle basis is 
limited to the lowest Landau level. The results are analyzed as a function of the total angular 
momentum of the system. Only at angular momenta corresponding to the filling factors 1, 1/3, 
1/5 etc. the system is fully polarized. The lowest energy states exhibit spin- waves, domains, and 
localization, depending on the angular momentum. Vortices exist only at excited polarized states. 
The high angular momentum limit shows localization of electrons and separation of the charge and 
spin excitations. 
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I. INTRODUCTION 

Semiconductor quantum dots have been a rich play- 
ground for both experimental and theoretical physics 1,2 . 
Nearly all methods of many-particle physics have been 
used to study the ground state properties and excitations 
in this well-defined problem of a few electrons trapped in 
a two-dimensional harmonic potential. The research has 
revealed the existence of a shell structure, Hund's first 
rule and many other properties related to atomic physics. 

In the presence of strong external magnetic fields the 
quantum dot is a finite-size realization of the quantum 
Hall liquid (QHL). In fact, it is one of the three ge- 
ometries (sometimes called the disc geometry) for per- 
forming exact many-particle calculations for QHLs. The 
other two geometries are the surface of a sphere, and a 
toru o 5 ' 6 : 7 . In a strong magnetic field, the electron system 
becomes polarized. Beyond magnetic field strengths cor- 
responding to filling factor v = 1 in QHL (or equivalently, 
to the maximum density droplet in a circular quantum 
dot) the spin does not usually play a role. Neverthe- 
less, several interesting phenomena, like edge reconstruc- 
tion^ and the formation of vortices&iiiSiiiii^ as well as 
magnetic excitations 1 ^*^ were observed. 

The Zeeman effect can be effectively diminished by 
choosing the material parameters such that the effec- 
tive Lande factor is zero. In so-called bilayer QHLs, a 
pseudo-spin gives the spin-like internal degree of free- 
dom (without the Zeeman effect). The spin degree of 
freedom drastically changes the 'simple' excitations of 
the QHL. The integer and the fractional quantum Hall 
systems at v — 1/3, 1/5 etc. stay polarized^ in their 
ground state and are thus called ferromagnetic, while 
at other fractions, like v = 2/3 and v = 2/5, the 
ground state is a spin singlet (S = 0). The QHL 
with spins has been analyzed with numerous theoretical 
methoda 5 ! 16 ! 17 ! 18 ! 19 : 20 ! 21 ! 22 ! 23 ! 24 ! 25 ! 26 . Many of these stud- 
ies are aimed at infinite QHLs, using periodic boundary 
conditions and (necessarily) restricting the numerical di- 
agonalization of the full many-body Hamiltonian to a few 



particles. 

The purpose of this paper is to study systematically 
the properties of many-particle spectra of simple sys- 
tems where a small number of electrons is confined in a 
two-dimensional harmonic potential. We solve the many- 
particle problem using the configuration interaction (CI) 
technique, the only approximation being that the single- 
particle basis is restricted to the lowest Landau level. 
Since we are mainly interested in large values of the total 
angular momentum (corresponding to small filling fac- 
tors) this approximation is suitable. The results show 
that many features of the spectrum are insensitive to the 
number of particles in the system and also to the specific 
form of the inter-electron interaction (remember that the 
Laughlin AnsataS 7 . for the fractional Hall effect does not 
contain any information about the Coulomb interaction) . 

In Section II, we will first describe the theoretical 
methods used. In Section III we discuss the results for 
the ground states and low energy excitations. 

II. THEORETICAL METHODS 

We assume a generic model of interacting electrons in 
a two-dimensional harmonic potential. The Hamiltonian 
is 

h 2 N N i N p 2 
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where N is the number of particles, m the electron mass, 
r = (x, y) a two-dimensional position vector, and ujq the 
oscillation frequency of the confining potential. Alterna- 
tively, we can use polar coordinates, x = 5i[r exp(z</>)], or 
complex coordinates z = x + iy. Note that we do not 
explicitely consider an external magnetic field. Instead, 
we solve the many-particle problem for a fixed angular 
momentum. In the case of a vanishing Zeeman splitting, 
the only effect of the magnetic field is to increase the 
orbital angular momentum of the system. 
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In addition, we also consider a contact interaction be- 
tween the electrons. Note that trivially, in this case only 
electrons with opposite spins interact: The Pauli exclu- 
sion principle forbids two electrons with the same spin to 
be at the same point. 

The Hamiltonian Eq. 0) is solved by using the config- 
uration interaction (CI) method, with the single-electron 
basis consisting of harmonic oscillator states in the lowest 
Landau level (LLL) 
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where Ag is a normalization factor and £ the single- 
particle angular momentum. In this basis, the diago- 
nal part of the many-particle Hamiltonian is constant, 
[L + N)huja, L being the total angular momentum. It 
is then sufficient to diagonalize the interaction part only. 
Therefore, the role of the strength of the confining poten- 
tial, ujq, is only to set the energy scale. In the LLL the 
structure of the many-particle spectrum and the many- 
particle wave function are completely independent of both 
luq and the strength of the electron- electron interaction 
(e 2 /47reo)- The restriction of the single-particle space to 
the LLL yields the exact result only in the weak inter- 
action limit. However, independent of the strength of 
the interaction, this approximation becomes more and 
more accurate when the total angular momentum of the 
system increases^. For the numerical calculation of the 
Coulomb matrix elements, see Stone et alm^. 

We diagonalize the Hamiltonian matrix for a fixed to- 
tal angular momentum L and a fixed z-component of the 
total spin, S z . We do not fix the total spin, but resolve 
its value from the eigenstates by taking the expectation 
value of the S" 2 operator. By selecting S z = for even 
number of electrons (S z = 1/2 for odd) the diagonaliza- 
tion of the Hamiltonian gives the energies and eigenstates 
for all possible values of the total spin. The dimension of 
the Fock space increases fast with N and L. For example, 
for N — 4, L = 30 the matrix dimension is 1234, but for 
N = 7, L = 42 it is already 43600. For large systems we 
make a further reduction of the Fock space by restrict- 
ing the maximum single-particle angular momentum of 
the space, £ < £ max - The numerical diagonalization was 
made using the Lanczos method 30 . It gives us the ground 
state and low-lying excited states. 

The symmetry of the system requires that the total 
electron densities and spin-densities are circularly sym- 
metric for all the states. In order to study the internal 
structure, we thus have to determine correlation func- 
tions. Here we use spin-dependent pair-correlation func- 
tions 



3TCT (r,r') = (*|n T (rK(r')|*), 



(3) 



where is the many-particle state in question, a de- 
notes spin up (t) or down (J,), and n-f is the spin-up- 
density operator. We define the total pair-correlation as 



g(r,r') = -(g n (r, r') + g n (v, r')) 



(4) 



The total angular momentum L can be related to the 
filling factor v via the relation 



N(N- 1) 
2L ' 



(5) 



For small particle numbers, this relation is strictly valid 
only for filling factors v = 1, 1/3 , 1/5, i.e. for 
those states of the fractional quantum Hall effect, which 
can be approximately described by the Laughlin wave 
functional. 

Within the LLL, the smallest possible angular mo- 
mentum is L2 — N(N/2 — l)/2 corresponding to a sin- 
gle Slater determinant where the spin-up and spin-down 
states are occupied from £ = to £ = N/2 — 1 (for 
even number of particles). It is natural to associate this 
state to filling factor v — 2, while Eq. (JjJJ would give 
v 2(1 + 1/(N+ 1)). Similarly, for fractional filling fac- 
tors 2/3 and 2/5 etc., Eq. (JSJ gives only an estimate of 
the corresponding angular momentum. 

We call the state with filling factor one the maximum 
density droplet (MDD)— . In this case, the angular mo- 
mentum is well defined, £mdd = N(N — l)/2. 

At large angular momenta the electrons crystallize in 
a rotating Wigner molecul e) 32 ' 33 . In this case the charge 
excitations can be described by classical vibrations of the 
molecule, separated from the spin excitations of the sys- 
tem. The equilibrium positions of classical electrons de- 
pend on the angular velocity oj r or angular momentum 
L = Iuj r of the Wigner molecule (/ is the moment of in- 
ertia I = Y] mr 2 ), and they can be solved by minimizing 
the classical energy 
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The eigenfrequencies of the vibrations can then be solved 
from the equations of motion of the rotating frame (by 
linearizing the equations around the equilibrium posi- 
tions of electrons )S. Quantization of the rotational and 
vibrational modes gives an estimate for the energy spec- 
trum 



Eqm = E cl 



1. 
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where lu^ are all the vibrational frequencies determined 
in the rotating frame and = 0, 1, 2, ■ • • , and the last 
term corresponds to the center of mass excitations. 

Once localized, the electrons form a system where the 
charge and spin excitations begin to separate. This is 
most clearly seen in one-dimensional systems^. In the 
case of quasi-one-dimensional quantum rings the whole 
many-particle spectrum can be quantitatively described 
with a model Hamiltonian 3 ^ 3 ^ 
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J^St-Sj , (8) 
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which consists of rigid rotations, vibrational modes of the 
localized electrons and of the Heisenberg model for the 
spin configuration. 



III. RESULTS 

A. Oscillating towards ferromagnetism, 2 > v > 1 

Figure 2] shows the total spin for the system of N = 
6 and N — 10 electrons as a function of the angular 
momentum. The spin increases, oscillating between zero 
and its maximum values, up to S — N/2, which is the 
spin of the maximum density droplet (MDD). Note that 
at this point the electron system is fully polarized (i.e., 
ferromagnetic). The lowest possible angular momentum 
for N = 10 (in the LLL) is L = 20. The corresponding 
state can be written as 
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where we denote the coordinates with the spin-down elec- 
trons as z. In the occupation number representation we 
write this 'double MDD' as 
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where the arrows show the occupied states of the LLL 
with increasing order of the single-particle angular mo- 
mentum. The wave function of the MDD, with S z = 
N/2, is exactly the Laughlin state 



MDD 
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which in the occupation number representation is 
I ^mdd) 
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The MDD has a spin degeneracy 2S + 1. For a given S z 
the state can be written as 

N/2+S z N/2-S z 

*mdd= Yl ( Zi ~ z i) II ( z k~zi) n(z m -z„)e~^ |z| 

i<j k<l m.n 

(11) 

It is important to note that for S z < N/2 the above wave 
function is a linear combination of several Slater deter- 
minants. For example, for S z = 0, these determinants 
are of the type 



TootoTTtoooooo 
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i.e. there is one electron at each angular momentum from 
to N — 1, but its spin can be up (f) or down (4). It 
is interesting to note that in the LLL approximation the 



exact wave function of the MDD is the (ferromagnetic) 
Laughlin state, although one would expect that much 
stronger correlations could be obtained by also allowing 
configurations with two opposite spins at the same angu- 
lar momentum. 

It is now easy to understand the maxima before the 
MDD, see Fig.Q] They consist of two maximum density 
droplets of different size, for example 



*(L=24)) 
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These states were studied in also in Ref« The states in 
between the maxima are more complicated and can no 
longer be described by a single determinant. For exam- 
ple, for L = 40, the most important configuration 
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has only 1.5 % weight in the exact wave function. 

Figure ^ shows that at small angular momenta (large 
filling factor) the system has a tendency to ferromag- 
netism. The states between the MDD's can be under- 
stood as spin waves of the ferromagnetic system, as ex- 
plained in the following section. The oscillations with 
increasing angular momentum, until the ferromagnetic 
MDD is reached, see Fig. ^ are independent of the par- 
ticle number. 
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FIG. 1: Total spin of the lowest energy state as a function of 
the angular momentum for N = 10 and N = 6 particles. 

The results for the LLL in the region 2 > v > 1 are 
'universal', i.e. independent of the electron- electron in- 
teraction. We repeated the calculations with a contact 
interaction. The resulting dependence of the total spin 
on the angular momentum was found to be the same as 
for the Coulomb interaction. Moreover, the simple states 
for the spin maxima are identical, and those in between 
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are also very similar. Table[I]shows overlaps of the many- 
particle states (K^couiombl^doita)! 2 ) for different angular 
momenta. The results demonstrate that for each angular 
momentum the lowest energy state is the same, but the 
order of the excited states can depend on the range of 
the interaction. 

TABLE I: Overlaps of the wave functions calculated with 
Coulomb and contact interaction for 10 electrons. L is the an- 
gular momentum and Si the spin of the state. The subscripts 
and 1 refer to the lowest energy state and the first excited 
state, respectively. For L = 44 and 45 the overlaps marked 
by * are those between the first excited state for Coulomb, 
and the second excited state for contact interaction. 
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For the MDD the overlap equals one, since in both 
cases the state is the simple Laughlin state of Eq. 
(consisting for N = 10 electrons of 256 Slater determi- 
nants with the same weight). To appreciate the good 
overlap between the more complicated states, one should 
note that, for example, for L = 40 the 1000 most impor- 
tant Slater determinants only contribute 97.0 % of the 
total wave function. 



B. Spin waves as excitations of the ferromagnetic 
state 

The completely polarized MDD is the ferromagnetic 
integer quantum Hall statei^2i2S. When the angular 
momentum is increased, the polarization decreases lin- 
early to zero at angular momentum L = £mdd + N/2. 
At larger angular momenta, the total spin of the lowest 
energy state is small, usually or 1, until one approaches 
the filling factor v — 1/3 where the system again becomes 
polarized, as seen in Fig. \I\ior six electrons {y = 1 cor- 
responds to L — 15 and v = 1/3 to L = 45). 

Figure [5] shows the oscillations in the total energy for 
10 electrons as a function of the angular momentum. The 
curves are obtained by subtracting a third-order polyno- 
mial fit from the total energy. The oscillating part of the 
total energy, here in the interval L = 20 to L = 45, 
is similar for Coulomb and contact interactions. The 
downwards cusps correspond to the spin maxima. As 
explained above, these states consist of two maximum 
density droplets with different size for spin-up and spin- 
down electrons. 

To understand the nature of the other states, we first 
study the excitation spectrum at the vicinity of the MDD. 
Figure shows the energy spectrum for 10 particles from 
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FIG. 2: Total energy as a function of the angular momentum 
for 10 electrons interacting with Coulomb interaction (upper 
panel) or delta function interaction (lower panel). Black dots 
indicate the states with local spin maxima (see Fig. and 
the squares the S = states. The energy is obtained by 
subtracting a third order polynomial fit from the total energy. 



L = 40 to L = 50. To illustrate the symmetry of the 
spectrum around the MDD (L = 45), we subtract a 
linear function from the total energy: Figure [3] shows 
AE = E- Hu (L + 1) - 19.86 + 0.162L, where the last 
two terms are obtained by fitting a linear function of the 
interaction energy (between L — 40 and L — 50). The 
spectrum is shown for the Coulomb interaction. In the 
case of the contact interaction, the energy differences be- 
tween the different spin states disappear after the MDD. 
The reason is that the MDD wave function, Eq. I|10|) , has 
zero amplitude whenever two electrons are at the same 
position and, consequently, the interaction energy of the 
contact interaction is zero. Beyond the MDD, the wave 
function can be constructed by multiplying the MDD 
state with a symmetric homogeneous polynomial 4 ^. 

Figure|21has a remarkable similarity with results in the 
study by Wojs and Quinn&Si, who considered a small 
number of electrons on the surface of a sphere. A similar 
spectrum is also obtained when using the geometry of a 
torusi. The structure of the spectrum seems then to be 
independent of the boundary conditions, though in each 
case the meaning of the angular momentum is different. 

The polarized state is the lowest-energy state only at 
the MDD (and at the higher Laughlin states). For other 
angular momenta, the spin is reduced. 

The spin of the ferromagnetic MDD can be reduced 
by spin-wave excitations. The operator that excites spin 
waves can be written aa 1 ^ 

N 

E+ = £ v^TT c+ +u Q t (12) 
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FIG. 3: Energy spectrum around the maximum density 
droplet of 10 electrons. Different symbols denote different 
spin: S — 5, black dots; 5 = 4, blue filled squares; S = 3, 
red open squares; S = 2, blue stars; 5 = 1, green x; S = 0, 
red +. The lowest energies for the fully polarized states are 
connected with a line. 



where c is the normal fermion annihilation operator. 
Note that operating with T,f to the MDD the angular 
momentum is increased by one and the total spin is de- 
creased by one. Similarly, we can form spin waves which 
lower the angular momentum: 



N 



c t-u c ei- 



(13) 



Moreover, assuming the spin waves as independent ex- 
citations, as in the theory of ferromagnetism^i, we can 
create two or more spin waves by operating successively 
with the above operators, for example, the state with 
two spin waves and angular momentum L — Lmdd + 2 
is (without normalization) 
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Table ITU shows overlaps between the exact result and 
that of the spin wave approximation. For the single spin 
wave the results agree with those of Oaknin et alM. For 
N = 10 electrons, the overlap decreases with the number 
of the spin waves, but is still more than 60 % for four spin 
waves. However, for the singlet state, which requires five 
spin waves, the overlap is less than 30 %. Clearly, this 
singlet state has a different character. This state (L — 40 
for N = 10) also shows a small kink in the total energy as 
a function of the angular momentum (see Fig. above) . 
We will return to this state in the next section. Oaknin 
et alw^ have shown that by projecting the simple spin 
wave state, Eq. <|12|) . to a state which is orthogonal to 
the center-of-mass excitations, the overlaps become even 
better. 



TABLE II: Overlap between the exact result and the spin 
waves states. N, L, and S are the number of electrons, angu- 
lar momentum and total spin, respectively, n is the number 
of spin waves. 
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FIG. 4: Occupation of the single-particle levels of the yrast 
state at different angular momenta for 18 electrons. Red line: 
L = 153 (MDD), green line: L = 154 blue line: L = 155, pink 
line: L = 156. 



Tabic UTI shows that when the number of electrons in- 
creases, the approximation of independent spin waves 
gets more accurate. For N — 18 electrons, even the three 
spin waves describe the exact result with 96 % accuracy. 
Figure Q] shows the single-particle occupancy of the ex- 
act states corresponding to 1, 2 and 3 spin waves for 18 
electrons. Clearly, the occupancy of the minority spin 
states (spin-down) increases linearly as a function of the 
single-particle angular momentum (for small angular mo- 
menta) in agreement with the spin wave suggestion, Eq. 
I|12p. Moreover, the increase of the spin waves increases 
the occupancy linearly. These observations are in qualita- 
tive agreement with the theory of Doretto et aZ.-- 6 , which 
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shows that the magnetic excitons have bosonic nature. 
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FIG. 5: Occupation of the single-particle levels of the lowest 
S — 8 state at different angular momenta for 18 electrons. 
Red line: L = 154 (MDD), green line: L = 155 blue line: 
L — 156, pink line: L — 157. (The red line is the same as the 
lower green line in Fig. [I] Note the different scale). 

Increasing the angular momentum from the MDD, the 
energy spectrum shows a large energy gap between low 
energy states and the rest of the states (including the 
fully polarized state, see Fig. We demonstrated pre- 
viously that the lowest of these states consist of spin- 
waves. The highest state below the energy gap has spin 
S = N/2 — 1 in the whole region from L = Lmdd + 1 
up to L = Lmdd + N/2 (i.e. from L = 46 to L = 50 for 
N = 10). Figure |3 shows that the minority spin concen- 
trates closer and closer to the surface when the angular 
momentum increases. The total spin density n-\ — still 
remains positive at the surface. This suggests that the 
spins along the surface are slightly tilted 4 ^. These states 
are often referred to as having non-collinear (or canted) 
spins 4 ^ 3 -. We should note that in our computations we 
have chosen S 2 and S z as good quantum numbers. Con- 
sequently, we do not know the values of S x and S y . 



C. Domain walls 

Large ferromagnetic systems may form domains. In 
normal magnetic materials the domains are caused by 
the long-range dipole interaction and exist already in 
the ground state. However, even without the dipo- 
lar interaction, domains can form as low energy exci- 
tations. This seems to be the case in quantum Hall 
liquid :] 5 - 15 - 20 - 23 - 44 - 45 - 46 

In finite small systems, the existence of domains is not 
as obvious, due to the small number of electrons and the 
finite thickness of the domain walls. Nevertheless, our re- 
sults show a clear separation of up and down spins at cer- 
tain angular momenta. Figures and {7\ show a set of pair 
correlations showing the tendency to separate the spin- 
up and spin-down electrons. The gradual development to 



a 'domain wall state' is indicated in Fig. which shows 
the correlations for a dot with N = 10 electrons from 
angular momenta Lmdd to Lmdd + N/2. The largest of 
these angular momenta (L = 50) shows clearly, that the 
up and down electrons prefer to be at the opposite sides 
of the dot. (If the reference electron, say spin-up, was at 
the center of the dot, the other spin- up electrons would 
stay closer to the origin, while the spin-down electrons 
are pushed further out). 



j \ , 



L = 46 
MDD 



L = 48 



L = 49 L = 50 



FIG. 6: Pair correlation functions beyond the MDD for 
N = 10 electrons. The upper panel shows the up-down corre- 
lation g-fi, and the lower panel the up-up correlation, g^. The 
number of spin-down electrons increases from one to five, and 
the number of spin-up electrons decreases from 9 to 5 when 
going from left to right. 
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f 



FIG. 7: Pair correlations for spin-singlet states after the 
MDD. The upper row shows the up-down correlation, gi,i, 
and the lower row the up-up correlation, g-f.T- The number 
of electrons N and the angular momentum L are indicated. 
Note that N = 10, L = 55 corresponds to the state N = 8, 
L = 36. 

The structure of two domains for L = Lmdd + N/2 is 
a general feature and does not depend on the number of 
electrons. Calculations were performed for N = 6, 8, 10, 
and 12 electrons. In all cases, the development of the pair 
correlation functions appear to be similar to that shown 
in Fig. [fl] above. 

Increasing the angular momentum further results in 
more complicated domain structures. For angular mo- 
mentum L = Lmdd + N, four domains seem to form as 
seen in Fig. Qfor N = 8, L = 36 and N = 10, L = 55. 
For even larger angular momenta it is not possible to re- 
solve possible domain structures at these small particle 
numbers. At angular momentum L = 44, corresponding 
to the filling factor v = 2/3, the pair correlation function 
can not be as easily interpreted. 

It is important to note that our computations do not 
contain the dipole-dipole interaction of the magnetic mo- 
ments of electrons. The origin of the domain formation 
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is thus different from that in normal ferromagnets. An 
infinite system would not have domains in the ground 
state. The fact that they appear in the lowest energy 
state at a fixed angular momentum does not imply that 
they represent a ground state. The antisymmetry of the 
total wave function restricts the allowed quantum states 
of a given total angular momentum. For example, the 
polarized state has lowest energy only for a few angu- 
lar momenta, while others must have some other internal 
structure. This becomes clear in the limit of localized 
electrons, to be discussed later. 



D. Vortices and edge reconstruction 

In a strong magnetic field, the Zeeman effect polar- 
izes the electrons and the system usually stays polarized 
beyond the maximum density droplet. The low-energy 
excitations of the MDD are then characterized by vor- 
tice a 8 ! 9 ! 11 ! 12 or edge reconstruction^^. As discussed ear- 
lier by Yang et ai4£, when the angular momentum in- 
creases from that of the MDD, holes will be formed in 
the otherwise filled Fermi sea below the single-particle 
angular momentum i = N — 1. Close to the MDD, the 
most important configuration of the many-particle state 
will have only one hole. For example at L = 32, N = 8 
we have 
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Increasing the angular momentum, the hole stepwise de- 
creases its (single-particle) angular momentum, "moving 
to the left" in the Fock state until it reaches I = 0. In 
the full CI calculation for polarized electrons, there are 
many other possible configurations with the same angu- 
lar momentum, but this single-hole configuration has the 
largest weight. This is clearly seen in Fig. |SJ where 
the thick red line shows the result for polarized electrons 
(5 = 4). 

The formation of a sing le hole in the MDDii can be 
understood in two ways. We can consider it as a single- 
particle excitation, where an electron is excited from 
I = 4 to I = 8. Alternatively, we can view it as a col- 
lective excitation, where four electrons from I — 4 • • • 7 
are each excited to the next angular momentum state. 
The latter interpretation is consistent with the vortex 
picture in boson systems^ where the Bose condensate cor- 
responds to the MDD of spinless fermions. 

The holes in the electron sea of the MDD can be iden- 
tified with vortices^. The electron density has a mini- 
mum at the site of the vortices, the wave function phase 
changes by 2-7T in going around the vortex center— and the 
electron current circulates around the vortessSi. 

We now return to the case without Zeeman splitting. 
The total spin of the lowest-energy states beyond the 
MDD is then generally small, and the fully polarized state 
is a rather high-lying excited state. Figure 03 shows the 
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FIG. 8: Single-particle occupancy of different energy states 
for N — 8, L = 32, 34, and 36. The occupancies of the 
lowest fully polarized state 5 = 4 are shown by thick red 
lines, the states below the first polarized state are shown by 
thin lines. The polarized results show a clear dip at a given 
single-particle state as a signature of a single vortex, which 
reaches origin at L = 36. 



single-particle occupancies for all the energy states below 
the polarized states, at angular momenta L = 32, 34 and 
36. A typical dip in the occupancy at one single-particle 
angular momentum is seen only in the fully polarized 
state. The low-spin states show a more uniform reduction 
of the occupancy and softening of the surface. Vortices 
seem to be low-energy excitations only for the polarized 
electron gas. 
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FIG. 9: Upper panel: Occupancy of the single-particle an- 
gular momenta of the lowest state for N = 8, L = 42 (red 
line), and for the lowest 5 = 3 (blue line) and 5 = 4 (green 
line) states. The lower panel shows the hole-hole correlation 
functions for the same states. 



The most clear two-vortex state in the case of eight 
electrons is at angular momentum L = 42 (for polarized 
electrons) 9 . In this case, the most important configura- 
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tion is 
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More generally, for n vortices the most important con- 
figuration is the one where the n holes are next to each 
other 4 i Fi gure El shows the single-particle occupancies 
of the exact diagonalization for N = 8, L = 42 for all 
the states up to the fully polarized state. Again, we see 
that only in the polarized state, shown by the thick red 
line, there clearly are two adjacent states (£ = 1 and 
£ = 2) where the occupancy is small compared to the 
other states. The only other state which has a rather 
similar occupancy distribution, is the one with total spin 
S = 3, i.e. the almost fully polarized state, shown by the 
thick blue line. The other spin states show only a weak 
minimum at the region of £ = 1 or 2, but larger occupan- 
cies from £ = 10 up, i.e. at the surface of the electron 
cloud. In fact, the general feature of the occupancy dis- 
tribution for all the states is, that there is a maximum at 
£ = and a broad maximum at larger £, the center being 
at about £ = 7. This is the first indication that at high 
(total) angular momenta the electrons begin to localize. 
For N = 8 electrons, eventually one electron will localize 
at the center, while the other seven form a ring around 
it. 

The excited states with 5 = 4 and 5 = 3 show oc- 
cupancy distributions consistent with existence of two 
vortices 9 . Other signatures of vortices are the phase of 
the wave function, and the pair correlation functions. 
The phase is not easy to compute in our formalism. Fur- 
thermore, as it is a function of 2N coordinates its inter- 
pretation is not straightforward^. Manninen et al£& have 
shown that by interpreting the vortices as holes, their lo- 
calization can be seen clearly in the hole-hole correlation 
function. Figure E] shows the hole-hole correlations for 
the lowest energy state 5 — and for the 5 = 3 and 
5 = 4 states. The correlations for the high spin states 
appear very similar. The reference vortex (or hole) is 
seen as a deep minimum. At the opposite side, there 
is a maximum corresponding to the second vortex. The 
ring around is caused by the holes outside the electron 
distribution. (In reality, the hole density is constant out- 
side the electron distribution, but since here the available 
single-particle space is limited, the hole density also goes 
to zero) . The hole- hole correlation function for the 5 = 
state does not even show a clear minimum at the site of 
the reference point (it is filled by the hole with opposite 
spin). The minimum seen in the case of 5 = is at 
the center of the dot and is caused by the fact that the 
electron density there is large and, consequently, the hole 
density is small. 

For large electron numbers (N > 20) the vortices of 
the polarized system can be seen clearly in the many- 
particle spectrum. The localization of vortices causes 
periodic oscillations to the spectrum^ in a similar fash- 
ion as the localization of the electrons when their num- 
ber is small 32 . Unfortunately, we are not yet able to do 



accurate computations for such large systems, when in- 
cluding the spin degree of freedom. However, the above 
results for N = 8 electrons strongly suggest that while 
the vortices are low-energy excitations for the polarized 
electrons, they are not the lowest excitations for electrons 
with spin. 

In the case of polarized electrons, the MDD breaks up 
by separating a ring of electrons^*. In the single-particle 
occupancy, the formation of this so-called Chamon-Wen 
edge^ is seen as a minimum in a similar way as the for- 
mation of vortices. Indeed, these two phenomena are 
intimately related. The minimum between the separated 
electron ring and the rest of the MDD is actually caused 
by one or two vortices^ ^ 9 -. In the circularly symmetric 
electron density, the ring of vortices is seen as a minimum. 
In the vortex-vortex (hole-hole) correlation, however, the 
vortices are localized along this ring. Again, at present 
we can not perform computations for large enough elec- 
tron numbers to study the possibility of Chamon-Wen 
edge reconstruction in the low-spin states. However, due 
its close relation to the vortex formation, we can safely 
predict that it will only be formed in the excited states 
with high total spin. 
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FIG. 10: Spectrum for N = 6 electrons for filling factors be- 
tween v = 1 and v = 1/3, as a function of angular momentum. 
Blue crosses: S — 0; blue squares: S = 1, red open squares: 
5 = 2, black bullets: 5 = 3. A third order polynomial fit- 
ted to the lowest energies has been subtracted from the total 
energy. 

The many-particle excitation spectrum shows a re- 
markable reflection symmetry around the angular mo- 
mentum L — N(N — 1), corresponding to filling factor 
v = 1/2. The composite fermion picture of Jain£2i£i 
reproduces this symmetry^ in a very natural way. Sub- 
tracting a smooth function of angular momentum, F(L), 
from the original spectrum removes the downward slope 
of the interaction energy, so that the small energy dif- 
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FIG. 11: Spectrum for 7 electrons v = 1 • • • 1/3. Blue crosses: 
5 = 1/2; blue squares: S 1 = 3/2; red open squares: S = 5/2; 
black dots: S = 7/2. A third order polynomial fitted to the 
lowest energies has been subtracted from the total energy. 
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FIG. 12: Spectrum for 8 electrons v = 1 • • • 1/3. Blue crosses: 
S — 0; blue squares: S = 1; red open squares: S = 2; red 
triangles: S = 3; black dots: 5 = 4. A third order polynomial 
fitted to the lowest energies has been subtracted from the total 
energy. 



ferences between the excited states can be more clearly 
illustrated, and the symmetry around the u = 1/2 states 
is mapped out in a very clear way. Figures 1101 1111 and 
IT21 show the spectra for 6, 7, and 8 electrons. In these 
cases, the function F(L) is a third order polynomial fit- 
ted in each case to the lowest energy states in the in- 
terval shown. At filling factor 1/2, the lowest state has 
total spin 0, 1/2 and 1, for N — 6, 7, and 8, respectively. 
At L = N(N — l)/2 and L = 3N(N - l)/2 (v = 1 and 
v = 1/3) the lowest state is ferromagnetic (has maximum 
spin). There are two other energetically favorable states, 
at L = N(N-l)/2+N 2 /±axi&L = 3N(N-l)/2-N 2 /4. 
Our interpretation is that the former of these corresponds 
to the filling factor v = 2/3 and the latter to v = 2/5 
(note that Eq. JSJ does not give exactly the above angular 



TABLE III: Overlap between the exact result and the 
Halperin-Haldane model, (\£/hjh | > 3/) | 2 , for different filling fac- 
tors v and electron numbers N 
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momenta for these special filling factors). The structure 
of the spectra around these special points is independent 
of the number of electrons (for an odd number of elec- 
trons, the lowest spin is naturally S = 1/2). 

The maximum spin state, S = N/2, is the lowest en- 
ergy state only at angular momenta corresponding to the 
filling factors v — 1 (MDD) and v = 1/3. It is interesting 
to observe that while the MDD has a large energy gap to 
the first excited state, the filling factor v = 1/3 does not. 
In fact, the energy gaps of the special singlet states at 
filling factors v = 2/3 and 2/5 have much larger energy 
gaps between the ground state and the lowest excited 
state. The large excitation gaps at these filling factors 
are well-known, see for example Ref. 7 . 



F. Halperin-Haldane wave functions 

The clear local minimum in the energy at angular mo- 
menta L = N (N - l)/2 + N 2 /A and L = 3N(N - l)/2 - 
N 2 /A suggests that there is a special way to build corre- 
lations in the wave function. The Ansatz wave function 
for the fractional quantum Hall effect 27 , Eq. (|ll|l , can be 
extended for some of the non-simple fractions as^iSi 
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where q is an odd integer and p a positive integer which 
can be even or odd. The angular momenta L — N(N — 
l)/2 + AT 2 /4 andL = 3N{N- l)/2-N 2 /A agree with the 
above wave function with q = 1, p = 2 and q = 3, p = 2, 
respectively. Table lTTTl shows the calculated overlaps. For 
the L = N(N — l)/2 + A^ 2 /4 state the overlap decreases 
rapidly with increasing N. Clearly, the exact quantum 
state can not be described by Eq. i|15|) for v = 2/3. In 
fact, in this case the wave function Eq. (|15f) is a mixture of 
the S = and S = 2 states. For L = 3N(N-l)/2- N 2 /4 
the agreement is much better, indicating that the above 
wave function is a good approximation to the ground 
state of v = 2/5. Also, it now has 5 = 0. 

It should be noted that for smaller filling factors than 
1/3 (or 1/5) the above simple analytic Ansatz does not 
accurately describe the true electron state. In particu- 
lar, it fails to describe the electron localization. In this 
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region, other trial wave functions, like for example the 
ones suggested by Yannouleas and Landman^, seem to 
be a better alternative. 

G. Symmetry of the spectrum around v = 1/3 
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FIG. 13: Spectrum for N = 7 electrons around angular mo- 
mentum L — 63 (y — 1/3). Blue crosses: S = 1/2; blue 
squares: S = 3/2; red open squares: S — 5/2; black bul- 
lets: S = 7/2. The upper panel shows the energy spectrum 
with a smooth part subtracted. The lower panel shows the 
order of different spin states in the low-energy band, below 
the large energy gap in the spectrum. Note the symmetry 
around L = 63. 

In figure above, we showed that after subtracting 
the smooth part, the resulting energy spectrum is sym- 
metric around the MDD (y = 1). It is thus interest- 
ing to see if similar symmetry exists also around the an- 
gular momentum L — 3N(N — l)/2, corresponding to 
v = 1/3. Figure ITol shows the excitation spectrum for 
7 electrons, now plotted such that the center is around 
L = 63. Again, a smooth function (in this case, a sec- 
ond order polynomial fitted to the lowest energies in the 
angular momentum interval shown) is subtracted from 
the original spectrum. The figure shows that in the im- 
mediate vicinity, from L = 58 to L — 68 the spectrum 
shows a qualitative mirror symmetry for the states in the 
low-energy band, below the large energy gap. The lower 
panel of Fig. 1131 shows the order of different spin states 
for the low energy band. We can see that, in addition 
to the symmetry near L — 63, there is a repeating pe- 
riod of six in the appearance of different spin states in 
the whole range show (ignoring their order). This is due 



to the emergence of the localization of the electrons in 
a six-fold ring with one electron at the center, as briefly 
discussed in the next section. 



H. Localization of electrons at high angular 
momentum, i/<l/3 

Increasing the angular momentum, the electron cloud 
in the harmonic confinement expands and the electrons 
begin to localize. For spinless electrons, it is possible to 
perform accurate CI calculations for much higher angu- 
lar momenta than if spin is included. For small electron 
numbers, the spectrum shows a characteristic periodicity 
as a function of the angular momentum2*22i22i2ii£§i£I. In 
these studies the number of electrons has been so small, 
that filling factor down to v = 1/9 have been reached. 
The results are then consistent with the expectation that 
in quantum Hall liquids the crystallization is expected^ 
to happen at filling factors smaller than v = 1/7. The 
periodicity in the spectrum as a function of angular mo- 
mentum is determined by the symmetry group of the 
Wigner molecule: Up to five electrons, the localized elec- 
trons form a single ring, the period being equal to the 
number of electrons. From N — 6 to 8 electrons, one 
of the electrons localizes at the center of the parabolic 
trap, while the rest form a single ring around : .t r '' u,i \ The 
length of the period of the oscillations in these cases is 
then N — 1. From figures ^2 and ^| we can see clearly 
that at high angular momenta, the fully polarized case 
has a low-energy state only at every 5th and 7th angular 
momentum, respectively. In the case of six electrons the 
period of 5 becomes clear only at higher angular momenta 
than those shown in Fig. ^3 due to the competition of 
two possible classical configurations of the electrons in 
the Wigner molecule 28 . 

For polarized (spinless) electrons the tendency for lo- 
calization is rather insensitive to the interparticle interac- 
tion. It appears for long-range Coulomb interactions, as 
well as for short range Gaussian interactions. Moreover, 
similar localization patterns occur for fermions as well as 
bosons? 2 . In all these cases, the effect of the localization 
is seen as periodic oscillations in the energy as a func- 
tion of the angular momentum. There is one exception, 
however. The exactly solvable model of harmonic inter- 
particle interaction does not show this periodicity 40,61 . 
The reason is the large degeneracy of the energy states 
and the fact that classically, all particles interacting by 
a repulsive harmonic potential are localized at the origin 
(otherwise, they would not be confined). 

Nikkarila and ManninenS have shown that at large an- 
gular momenta the whole many-particle spectrum of ro- 
tating polarized electrons can be quantitatively described 
by quantizing the classically determined vibration fre- 
quencies of the Wigner molecule. This suggests that the 
charge excitations (vibrational modes) and the spin ex- 
citations will separate in a similar fashion than in one- 
dimensional quantum rings^*^, in which case the system 
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FIG. 14: Excitation energies as a function of the angular mo- 
mentum. Points are results of the exact diagonalization and 
the lines from the model of classical vibrations. Blue points: 
S — 2, open circles: 5 = 1, green squares: 5 = 0. 



can be described by the Hamiltonian of Eq. (JHJ) . 

We determined the vibrational modes of classical elec- 
trons in a 2D harmonic confinement using a rotating 
frame, as in Ref. 33, and compared the results with those 
obtained from the full quantum-mechanical CI calcula- 
tion. 

For N = 4 electrons, the result is shown in Fig. ITTI 
The figure shows the excitation energies for each angu- 
lar momentum, i.e. the total energy of the excited state 
minus the lowest energy for the same angular momen- 
tum. The quantum-mechanical results approach the clas- 
sically determined energies when the angular momentum 
increases. Different spin states of the QM calculation are 
marked by different symbols. We can see clearly that 
the lowest energy state and each vibrational state has a 
periodic pattern of spin states. 
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FIG. 15: Total spin of the lowest energy state as a function 
of the total angular momentum for TV = 4 and N = 5. 

The periodic appearance of the different spin states at 



each vibrational mode in Fig. 1141 gives further support 
to the separation of the spin and charge excitations. As- 
suming the model Hamiltonian of Eq. (JHJ) one can use 
group theory 35 - 62-63 to analyze which spin state corre- 
sponds to which angular momentum. The Heisenberg 
coupling constant J determines the energy splitting of 
different spin states, and its sign determines if the elec- 
tron system is ferromagnetic or anti-ferromagnetic. In- 
deed, the obtained spin states agree well with the pre- 
diction of the model Hamiltonian. However, the simple 
model can not explain quantitatively the energy differ- 
ences between the different spin states. For example, for 
angular momenta 18, 22, 24, etc., the purely rotational 
band has spin states S = 2 and 5 = 0. If the Heisen- 
berg coupling would be ferromagnetic, the S = 2 state 
would have lower energy than the S = state, while for 
anti-ferromagnetic coupling, this would be the other way 
around. 

Figure ^] shows the total spin of the lowest energy 
states for 4 and 5 electrons as a function of the angular 
momentum. Clearly, we can see that neither of these 
cases show a simple period of N, as it would be expected 
from a simple Heisenberg model of localized electrons. 

In the case of one electron per lattice site the Hub- 
bard model approaches to the anti-ferromagnetic Heisen- 
berg models. Since here the electrons are localized in a 
Wigner lattice, one should expect that the same model 
works equally well (although now the electrons are local- 
ized not by external potential, but by the mean field cre- 
ated by the other electrons). For quasi-one-dimensional 
quantum rings this indeed was shown to be the cas o 35 ' 36 , 
as well as for electrons localized in the corners of trian- 
gular or square quantum dotsS 5 .. However, the result of 
Fig. EDshows clearly that this does not hold for electrons 
localized in a two-dimensional harmonic confinement. 
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FIG. 16: Energy spectrum for four electrons in a quantum 
dot (£min = 0, £ m ax = 11) compared to that of a quasi-one- 
dimensional quantum ring (£ m i n = 3, ^ max = 11) determined 
by restricting the basis as explained in the text. Blue points: 
5 = 2, open circles: 5=1, green squares: 5 = 0. 



We can mimic a quantum ring by restricting the single- 
particle basis from low and high angular momenta, i.e. 
^min < i < 4nax- Figure El shows part of the energy 
spectrum computed with this restricted basis compared 
to that of the full calculation. The result for the 'ring' 
shows an anti-ferromagnetic spin arrangement in agree- 
ment with the Heisenberg model^, while for the dot the 
result is different. Note, however, that in both cases we 
observe the separation of the lowest, purely rotational 
band, with the same sequence of spin states at each an- 
gular momentum. 




L ■ 34. 1st L - 34. 2nd L • 38, 1st L • 36. 2nd 

FIG. 17: Pair correlation functions for N = 4 at angular 
momenta L = 34 and L — 38 for the two lowest energy states. 
The upper row shows the up-up correlation, g^f, and the 
lower row the up-down correlation, fff,4' The position of the 
reference electron is shown as a cross. 

Figure IT7I shows examples of the pair correlation func- 
tions, Eq. (J3J), for four electrons in the purely rotational 
states. We can see that for the fully polarized states, i.e. 
the 1st state for L = 34 and 2nd state for L — 38, the up- 
up and up-down pair correlations are similar due to the 
symmetry of the ferromagnetic state (we show the S z = 
result), while for spin singlet states [S = 0) the pair cor- 
relations show a tendency for antiferromagnetism. It is 
important to note also that the two states for L = 34 
and L = 38 are similar, but of opposite order in energy. 
However, the energy differences are extremely small as 
seen in Fig. El 

For polarized electrons it is possible to perform an ex- 
act diagonalization for large angular momenta where the 
localization of electrons is clearly seen^. The spin degree 
of freedom, however, increases the available Fock space 
drastically. For example, for N = 7 or 8 electrons we are 
limited to study only the region up to v w 1/3. How- 
ever, we can see that the localization begins already in 
this region. Figure 1181 shows the pair correlation func- 
tion for 7 electrons at L = 51 corresponding to the filling 
factor v = 2/5. The total correlation function, Eq. 
shown in the upper row of Fig. 1181 shows clearly that 
the electrons begin to localize in a geometry where one 
electron is at the center, and the remaining six electrons 
form a hexagon around it. The lowest-energy state at 
this angular momentum is well separated from the ex- 
cited states and can be approximatively described by the 
Halperin-Haldane wave function as mentioned above. It 
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FIG. 18: Pair correlation functions for N = 7 at angular 
momenta L = 51 for the lowest energy state, and for the 4th, 
5th, 8th and 9th excites states. The uppermost row shows 
the total correlation function, g-f^ + gi,i, the center row the 
up-up correlation gt,t, and the the lowest row the up-down 
correlation, <?t,j.. The position of the reference electron is 
shown as a cross. 



is interesting to note, that in this state the electrons are 
not as well localized as in the excited states (4th, 5th 
and 7th state). The 8th state is above the large energy 
gap seen in Figs. 1111 and 1131 In this state, the electrons 
are not as clearly localized as in the lower-lying states. 
This is in agreement with the beginning separation of 
the spin excitations from the charge excitations: Below 
the large energy gap in Fig. El au the excitations are 
spin excitations, which do not markedly change the pair 
correlation function, while above the gap, the excitations 
include a charge excitation which necessarily changes the 
pair correlation. We checked the pair correlations for all 
the states from L — 51 to L = 75. In all cases, the states 
below the energy gap are similar, showing the localiza- 
tion of electrons which gradually becomes stronger when 
the angular momentum increases^. 

Figure ITH1 shows also the up-up and up-down pair cor- 
relations. These support the view that the low-energy ex- 
citations are mainly spin excitations. The three 5=1/2 
states below the energy gap (i.e. the 1st, 4th and 5th 
state) show different up-up and up-down correlations, 
while the total correlation is nearly similar. In the case of 
the lowest energy state, the up-up and up-down correla- 
tions are in agreement with the Halperin-Haldane state: 
The repulsion between opposite spins is larger than that 
between the same spins. 

The pair correlation functions shown in Fig. 1 181 do not 
show much similarity to those we obtained in the anti- 
ferromagnetic Heisenberg model. 

Note, that in the fully polarized case, S = 7/2, the up- 
up and up-down correlations, calculated for S z = 1/2, 
are identical, as expected. 
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IV. CONCLUSIONS 

We have studied rotational states of interacting elec- 
trons in a two-dimensional quantum dot confined by a 
circular harmonic potential. The many-particle Hamilto- 
nian was solved exactly by numerical diagonalization in 
a basis restricted to the lowest Landau level. The results 
were analyzed in terms of the total angular momentum 
and spin. Magnetic fields were not considered explicitely, 
and consequently, no Zeeman splitting was introduced. 

The smallest possible angular momentum in the LLL 
is a 'double maximum density droplet' with L = N(N — 
2)/4 corresponding to the filling factor v = 2. Increasing 
the angular momentum the total spin of the system in- 
creases oscillating between zero and the next maximum 
until it reaches the fully polarized ferromagnetic state of 
the maximum density droplet. This behavior was found 
to be the same for the long range Coulomb interaction 
and for a contact interaction. 

Beyond the MDD, the lowest energy states are spin- 
waves of the ferromagnetic state in agreement with pre- 
vious calculations with different periodic boundary con- 
ditions. Increasing the angular momentum further, the 
lowest-energy states seem to consist of ferromagnetic do- 
mains, while the total spin of the system is zero. 

For low total spin we did not find any vortices as low 
energy excitations, although for the polarized case they 
appear as the lowest energy states at certain angular mo- 
menta. 

The excitation spectrum shows a clear reflection sym- 
metry around angular momentum L = N(N — 1) corre- 
sponding to a filling factor v = 1/2. In the neighbor- 
hood of angular momenta L — N(N — l)/2 + N 2 /4 and 
L = 3N(n — l)/2 — 7i 2 /4, corresponding to filling factors 



2/3 and 2/5 respectively, the spectra look very similar. 
Both states have S = and a large energy gap to the first 
excited state. In fact, these energy gaps are the largest 
in the whole region from v = 1 to v = 1/3, apart from 
those in the immediate vicinity of filling factor v = 1. 
The lowest-energy state with filling factor v = 2/5 can be 
rather accurately approximated by the Halperin-Haldane 
generalization of the Laughlin wave function, at least for 
small numbers of electrons. 

In addition, the spectra are symmetric around angular 
momenta L = N(N - l)/2 and L = 3N(N - 1)2 cor- 
responding to filling factors 1 and 1/3, respectively. Es- 
pecially, the states around v = 1/3 seem to have similar 
spin wave excitations as the states around filling factor 
one. Very recently, Dethlefsen et al~ have studied the 
details of the excitation in the region of v « 1/3 with 
comparison to experiments. 

At large angular momenta, the electrons start to local- 
ize. In the case of four electrons, accurate results could 
be obtained up to angular momenta L < 46, correspond- 
ing to a filling factor smaller than v = 1/7. In this case 
the results showed a clear separation of charge-like and 
spin- like excitations. The charge excitations could be 
quantitatively explained by quantization of the classical 
vibrational modes of the localized electrons. However, 
contrary what we expected (for localized electrons) the 
spin-excitations could not be explained within a simple 
Heisenberg model. 
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